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Abstract 

We study the Lepton number violating like-sign dilepton processes hih 2 1*1* jjX and 
hih,2 — ► I I W X, mediated by heavy GeV scale Majorana neutrinos. We focus on the reso- 
nantly enhanced contributions with a nearly on-mass-shell Majorana neutrino in the s-channel. 
We study the constraints on like-sign dilepton production at the Tevatron and the LHC on the 
basis of the existing experimental limits on the masses of heavy neutrinos and their mixings U a M 
with a = v^v^v T . Special attention is paid to the constraints from neutrinoless double beta 
decay. We note that searches for like-sign e^e^ events at Tevatron and LHC may shed light on 
CP-violation in neutrino sector. We also discuss the conditions under which it is possible to extract 
individual constraints on the mixing matrix elements in a model independent way. 
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I. INTRODUCTION 



At present there are no doubts that neutrinos are massive particles mixing with each 
other. Moreover, according to the neutrino oscillation experiments, their masses are ex- 
tremely small while mixing is nearly maximal. In this respect neutrinos drastically defer 
from all other known particles, and this difference represents one of the pressing problems 
for theory. The famous see-saw mechanism showed the way to possible solutions of this 
problem via introduction of very heavy Majorana particles mixed with ordinary neutrinos, 
bringing in the observable mass hierarchy between the light neutrinos and other fermions. 
These heavy particles could be heavy Majorana neutrinos, as in the original formulation of 
see-saw, or other new particles such as neutralinos, heavy Majorana particles mixing with 
neutrinos in the framework of SUSY models with R-parity violation. The new heavy parti- 
cles may lead to observable effects beyond the light neutrino sector, since they could manifest 
themselves indirectly via their virtual contribution to processes involving ordinary particles. 
If the masses of some of these new particles are not extremely large, lying within the reach 
of the running of forthcoming experiments, they can be searched for directly among the 
products of colliding particles. 

Here we consider an extended see-saw scenario, including n species of SM singlet right- 
handed neutrinos v' R - = (v' R1 , ■ ■■i / Rn ), besides the three left-handed weak doublet neutrinos 
u Li = (y'hei v Liii u Lt)- The general mass term for this set of fields can be written as 

-h?MMv* + h.c. = ~KK C ) (" L T [ <C \ +hx - (1) 

2 2 \M T D M R )\v' R ) 

3 n 

= - o (2 rn v .T^Ui + mvrfvj) + h.c. (2) 

Z i=l j=l 

Here Ml, Mr are 3x3 and n x n symmetric Majorana mass matrices, and M D is a 3 x n 
Dirac type matrix. Rotating the neutrino mass matrix to the diagonal form by a unitary 
transformation 

U T M M U = Diag{m I/ i, • • • , m, 3+n } (3) 

we end up with 3 + n Majorana neutrinos with masses m Vl , • • • , m V3+n . In this scenario the 
light neutrinos have the mass scale M 2 d /Mr, and as a result there should be also heavy 
Majorana neutrinos (N) with mass scale Mr. Due to the fact that Majorana neutrinos are 
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also their anti-neutrinos, they can produce lepton number (L) violation by two units. One 
of the decisive processes for probing the Majorana nature of neutrinos is neutrinoless double 



Majorana 



beta (0^/5/3) decay which has been extensively studied in the literature 
neutrinos may also resonantly contribute to meson [3J] and r decays |4|. 

Another potential process to look for Majorana neutrinos is like-sign dilepton production 
h\h 2 — > l^l'^W^X at hadron colliders {|], si- m this paper we give a detailed analysis 
of like-sign dilepton production processes, evaluating the cross section of hih 2 — > l^l'^jjX 
and h x h 2 -> iH^W^X. 

The existing very stringent constraint from 0i>(3(3 experiments on the inverse effective 
mass t/'gjy/Mjv of a heavy Majorana neutrino N is frequently treated as the limit, which 
makes unrealistic the experimental observation of the above processes with e ± e ± pair in the 
final state. However, we note that CP-violating Majorana phases, if present in the neutrino 
mixing matrix elements U ea , are able to significantly soften or even completely evade the 
Ot//3/3-constraints on the e ± e ± -process production. We point out that this can be true even 
in the presence of only one heavy Majorana neutrino, due to the interference of heavy-light 
neutrino contribution to 0i//?/?-decay. 

The paper is organized as follows. In the next section we present our approach to calcula- 
tion of the cross sections of the hih 2 — > I I jjX and h\h 2 — > l^l'^W^X processes. In Sec. 
Ill we calculate the decay width of the heavy Majorana neutrino entering in the above cross 
sections. Sec. IV is devoted to a discussion of the existing limits on the parameters of heavy 
Majorana neutrino and their impact on the prospects for searches for like-sign dileptons at 
Tevatron and LHC Here we comment on implications of the 0^/3/?-decay constraints and 
CP-violation in the neutrino sector and discuss the possibility of extracting the heavy-light 
neutrino mixing matrix elements U e N, U^n, UtN in a model independent way. 

II. THEORETICAL FRAMEWORK 

The processes we study are like-sign dilepton inclusive production in high energy pp or 
pp collisions 

h x (Pi) + h 2 (P a ) - Z±(Z0 + l f± {l 2 ) + jj{W*) + X, (4) 
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(a) 



FIG. 1: (a) resonant and (b) non-resonsnt Majorana neutrino contributions to 
hi{P\) + ^(-fb) - ► l + (h) + l' + (h) + X process. In the diagrams we indicate the momenta 
of the corresponding particles. 

where jj and X denotes two quark jets and undetected hadronic states. The two leptons 
/ and V can have the same or different lepton flavors. In these processes the total lepton 
number L is violated in two units AL = 2. If one assumes the existence of a Majorana 
neutrino N, the process can be realized through the diagrams shown in Fig. HJ We focus 
on the resonantly enhanced diagram Fig. [H(a) with the nearly on mass-shell neutrino in the 
kinematical region of the studied processes. 

In the diagram Fig. QJa) there is a W boson in the final state, which decays to hadron or 
lepton states. Here we will consider the situation in which W~ decays to two quark jets, as 
shown in Fig. [T^,. According to the factorization theorem, we can express the cross section 
of the process as 



here q a /h\{ x a) an d qb/h 2 ( x b) denote the densities of partons inside each hadron, and 




.q a qb^l + l' + qkqi 



+ (x a ^ x b )}, (5) 



d a qaqb^i + i' + qkqi j s ^ ne differential cross section of the partonic subprocess 



q a (k a ) + q b (h) -> (h) + I' + W(P W ) -> 1(h) + q h ( qi ) + g z (g 2 ) 



(6) 



which has the form 



qkqi _ 



M q a q b ^ + l' + q k m 2 ^f5\k a + k b - h ~ h ~ ?1 ~ ft) 



d 3 li d 3 l 2 d?qi d 3 q 2 



(7) 



x 



(2tt) 3 2£i (2tt) 3 2£ 2 (27r) 3 2E gi (2vr) 3 2£, 
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The next step is to calculate the amplitude squared \M.\ 2 , which can be written as 



MP - leajMi, x ^^'^j^i^^l? + (*. - *»), (8) 



((i - Ml? + r° w Ml)((Pl - Ml? + Tl.Ml) 
with the tensors having the following forms: 



= 16(A£A£ + KK - gTK ■ hMqi + rftf ~ gTli ' ft) 



E™>nUinUi'n -n \ u an cn \ n i \ 

? [(t-^-mM ''^ 7 ( ' 2) - ((l " W 



X 



E 



t^nUinUvn 



v(h)Y^P R v c (l 2 ) - (Z x <-> l 2 ) 



N [{k - h) 2 - m 2 NS 
The contraction of the above two tensors yields 



t 



(9) 



(10) 



inJ>a p = 128 V FTy-^TT^ 19 ((ii ■ foXfe ■ <?2)(&a ■ Qi) + (k ■ h)(k ■ gi)(fc a • 92) 
+ (Zi • k a )(h ■ q2)(h ■ gi) + (Zi • &b) (i 2 • q2){h ■ q 2 )) + interference term. (11) 



We have checked that the interference term gives a very small contribution compared to the 
first term in ( ITT]) , and therefore we ignore it in the following calculations. The 4-body phase 
space in Eq. (J7|) can be rewritten as 

d 3 l\ d 3 l 2 d 3 Pw 



$4 = S\k a + k b -h-l 2 -P, 



(2n)*2E h {2Tif2E l2 (2tt) 3 2E w 



x 5 4 (P W -qt- q 2 ) 



d 3 qi d 3 q 2 



WdP^, 



(12) 



{2i,f2E qi {2nf2E q2 

where Pw = qi + ?2 denotes the momentum of the final virtual W boson. In the cm. frame 
of qi and q 2 , the two body phase space can be expressed as 

d 3 qi d 3 q 2 



$ (Pw - <?i - ft) 



(2n) 3 2E qi (2ir) 3 2E 



X^ 2 (P w ,m 2 qi ,m 2 q2 ) ^ 
8(2^ P w 



(13) 



where dQ = dcos9\d(f)i is the solid angle of qi, and the form of \(P w ,m qi ,m 2 2 ) is given in 
Eq. flA5j) . In the case P^ 3> m 2 and m 2 , X 1 ^ 2 (P W: m 2 1 ,m 2 2 ) ~ P w . After integrating over 
cZfi, Eq. (ED turns to 
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3 S Kr^i'^ tft-w^+fa-*-) 

+ (k-k a ){2E l2 E kb + l 2 -k b )), 



(14) 



and the four-body phase space in Eq. (TT21 is reduced to a three-body phase space. 



As shown in Fig. [21 it is convenient to use the following Lorentz invariants in order to 
express the kinematical variables: 



(h + Pw) 2 , 



Sb = {h + Pw) 2 , t a = (k a - l 2 ) 2 , t b = (k b - h 



(15) 



With the above invariant variables, we can apply the phase space transformation given in 
Ref. {9J to rewrite the cross section given in Eq. (jSJ) as 

a = a dx a / dx b / ds a / dt a / & ds b j dt b dP w 

J J Js~ Jt~ Js~ Jt~ J 

x / a -M.{x a , Xb, s a , Sb,t a ,tb, P\y)j 

whith a = 2G 4 F M w /(3(2n) 6 ), and 

\V ud \ 2 u(x a )d(x b ) + \V us \ 2 u(x a )s(x b ) + (x a <-► x b ) 



(16) 



hA. (yX a , x b , s a , s b , t a , t b 



{{s - m&Y + t* w m&)((p& - M^y + r^) 
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(17) 



E 



N s b m N 



M i(s, s a , Sb, t a , t b , P w ) + interference term, 



where 



Mi(s,s a ,s b ,t a ,t b ,Pw) = {s - s b + t a -ml)((s b - P"^ -m 2 l2 )(s a + t a -t b -m 2 h ) 

+ P^(s - s a + t b - m 2 b )) - (t a -m 2 a - m 2 h )((s b - Pw - m 2 2 ) 



x (s b -t a + t b - mfj - P w (t b -m 2 b - mf 2 )), 



(18) 



The phase space can be constrained to the physical region by the following inequalities for 
the Gram determinants 



A 3 (Pl,P2,P3) > 0, A 4 (pi,p2,P3,P4) < 0, 



(19) 



FIG. 2: Definition of the invariant variables. 



where {pi} denotes any subset of momenta {k a , kb,h,h, Pw}- The resulting integration 
limits s~ /b , s+ /b ,t~ /b and £+ /6 are given in Appendix. 

If there exists one Majorana neutrino with the mass is in the region s b < m 2 N < s b , 
the integrand shown in Eq. ffTB]) has a pole at s b = m 2 N . Introducing a decay width for 
N through the substitution m N — > m N — (i/2)T N , the Majorana neutrino propagator can 
be written as l/(s — m 2 N + 277^1^). As we will show in the next section, -C for m^v 
in range we are studying in the present paper. Therefore, the cross section can be written 
as 



a hlh ' 2 ^ l+l ' +x (m N ) ps air \ dx a \ dx b \ ds a \ dt a j dt b \ P^ 



9(-A 4 ) 
s 2 ^Kl 



\V U d\ U(X a )d(Xb) + \V US \ U(Xg)s(X b ) + (Xg <~> X b ) 

X {{s-M 2 w ) 2 + T 2 w M 2 w ){{P^-M 2 w ) 2 + T 2 w M 2 w ) 

x mNlU ' NUl ' Nl2 M 1 (s,s a , m %,t a ,t b ,P^ (20) 

J- N 



To arrive at the above equation we have used the identity: 

f(s b ) _f(m N ) ^ 2 



ds b - Tv?~, 2~F2~ = n Fi — > s b < m JV < s 6 an d « m N . (21) 



(s 6 - m? N ) 2 + m^r^ m^r 



N 



Therefore, with in the denominator, the cross section can be strongly enhanced through 
the resonant production of a Majorana neutrino. If there are several Majorana neutrinos in 
the range < m 2 N < s b , there should be a sum over Ni in Eq. f l20l . 

In the case > Mw, one can detect the W boson together with the lepton pair 
in the final state. Therefore we also give the formula to calculate the cross section of the 
process h x + h 2 -> 1+ + l'+ + W~ + X, 

a h ^ l+l ' +w - x (m N ) « a'n [ dx a [ dx b f ds a f dt a f dt b ^t^ 



x 



t~ Jt b s 2 x/-A 4 
V ud \ 2 u(x 1 )d(x 2 ) + \V us \ 2 u(xi)s(x 2 ) + (xi <-> x 2 ) 



x ^ L jM' 1 (s,s ,m^,t a ,t 6 ), (22) 

J- N 



where 



M' x (s,s a ,s h ,ta,t b ) = (s a + t a -t b -m 2 J(s - s b + t a -ml)(s b - -mf 2 )/Mw 

- (t a -m 2 a - ml)(s b -tg + t b - m 2 3 )(s b - - m 2 h )/M^ 

+ (ml + ml - t b ){m 2 a + m\ - t a ) 

+ (s- s a + t b -ml)(s- s b + t a -m 2 a ), (23) 



and a' = 2 v / 2G^M^/(6(27r) 4 ). 



III. NEUTRINO DECAY WIDTH 

As shown in Eq. (12Uj) and (1221) . the cross sections of the process hih 2 — > Z + Z /+ X, and 
hih,2 — > / ± / /± W /:F A mediated by the resonant production of Majorana neutrino depend on 
the neutrino decay width T^. For m N <C M^/ the decay width Y N receives contributions 
from leptonic and semi-leptonic channels 

N^hhu, Vi{vt)U, l^qxqi, Vi(v£)qq. (24) 

We approximate the semi-leptonic decays by inclusive quark-antiquark pair production. In 
such an approach {4] based on the Bloom-Gilman duality [lo| the total decay width of the 
heavy Neutrino N to hadrons is reproduced in average with an accuracy sufficient for our 
estimates. One of the advantages of this approach is that it does not require knowledge 
of meson masses M M and decay constants F M necessary for the calculation of N — > IM. 
partial widths which then are summed up in order to derive the total decay width of heavy 
neutrino N in the channel-by-channel approach . For heavy mesons these parameters are 
purely known and may introduce a significant uncertainty. 

The leading order decay rates for the channels listed in (j24"|) can be found for the leptonic 
and semileptonic decays in Refs. Q, fl 3 and Q respectively. We summarize the corre- 
sponding formulas neglecting lepton and quark masses. In this approximation the partial 
decay widths of Majorana neutrino in the region < My/ are 

T(N - hl 2 u k ) = \U hN \ 2 ^^m 5 N F w (m N ) = \U hN \ 2 T l {, (25) 
T(N - v h l 2 h) = \U hN \ 2 ^m%F z (m N ) x (26) 

x [{g l L ? + (g l R ) 2 + qWr + <U(ak + = \u hN \ 2 ^, 

£ V{N^v h V l2 v l2 ) = \U hN \ 2 ^m%F z {m N ) = \U hN \ 2 r 3 », (27) 



l2=e,fi,T 



T (N - l^ud) = \U hN \ 2 (|Kd| 2 + \V US \ 2 + \V cs \ 2 )-^m%F w (m N ) = (28) 



64vr 3 



I t t 1 2j^lud 
\ U hN\ 1 , 



G 2 

T(N^u iq q) = \U hN \ 2 ^m 5 N F z {m N ) x (29) 

x M) 2 + (g R ) 2 + g q L g R ] = \u hN \ 2 r^, 

8 



In Eq. fl28l) we neglected the small charged current cb contribution. Here g l L = 
-1/2 + sin 2 6 W , g l R = sin 2 9 W , g u L = 1/2 - (2/3) sin 2 W , g\ = -(2/3) sin 2 9 W , g% = -1/2 + 
(1/3) sin 2 Q w , g\ = (1/3) sin 2 Q w are the SM neutral current lepton and quark couplings 
where I = e, /x, r and u — u,c,t; d — d, s, b. 
In Eqs. (j25) - (E9|l the function 

F B (m N )=6M B J Q d Sl J Q ds 2 <((s2 _ M|)2 + r|M|) • (30) 

takes into account a propagator effect of W and Z-bozon exchange and introduce a significant 
correction of a factor ~ 3 for tjin ~ Mw, Mz- We denote B = W, Z. Here Mw = 80.4GeV, 
T w = 2.14GeV and M z = 91.2GeV, T z = 2.5GeV are the mass and full decay width of W 
and Z-bosones respectively. 

Therefore, in the case rriN < Mw, the total decay rate is 



T N (m N < m w ) ps ( J! \Uin\ 



|2 



2T[ + 2 r lud + r% + t 3u + vvq 

q=u,d,s,c,b 



• (31) 



The factor 2 in the first two terms is due to Majorana nature of heavy neutrino N which 
can decay into two non-equivalent charge conjugate final states. For m N > M w the total 
decay of heavy neutrino is determined by the decay channels: 

N -> l T W T , N — > uZ, N — > vH°. (32) 



The corresponding partial widths of Majorana neutrino are given by [6|, 111! ]: 

r(iv-^) = |^<(i + ^)(i4)^K-%)Har w , 

8V27T m z N m z N 

(33) 

r(tf - *>(if)Z°) = |f/ w | 2 ^m5(l + ^)(l-% 2 ^-M z ) = |^| 2 r^) ; (34) 

8V27T J m^r 

r(iV - ^)#°) = \U lN \ 2 -^rm%l - ^r) 2 d{m N - M H ) = \U lN \ 2 T^ H \ (35) 

8V27T J m z N 

In the same way we can write the total decay width of a Majorana neutrino at mjy > Mw 

as 

T N (m N > M w ) ~ ( E l^l 2 ) ( 2rM/ + rZ + TH ) (36) 

\/=e,/i,T / 

The above formulas we use in our analysis of the like-sign dilepton production. Numeri- 
cally r^v <?C TTiN holds for lGeV < < ITeV and, therefore, Eqs. (|20|) and (1221) are good 
approximations for the cross sections of the process h 1 h 2 — > / + / /+ jj(W / )X. 



IV. LIMITS ON HEAVY MAJORANA NEUTRINO SECTOR 



The remaining ingredient, which determines in the most crucial way the event rate of the 
like-sign dilepton production, is the heavy Majorana neutrino masses and their mixing with 
v e , v T . In the literature there are various limits on M N and U a N for a = e, /i, r, extracted 



from direct experimental searches [12j for these particles in a wide region of masses, from the 
non-observation of lepton number violating decays, as well as from precision measurements of 
certain observable quantities [3]. In the present context we are interested in these limits for 
the mass range > lGeV. From the global fit of the electroweak precision measurements, 
including LEP data, it was found [3]: 

\U eN \ 2 < 0.0052, |[7^| 2 < 0.0001, \U tN \ 2 < 0.01 (37) 

The LEP data on direct searches of heavy leptons [14J] shows that 

\U tN \ 2 < 10" 4 - 10" 5 for 3GeV < M N < 80GeV (38) 

where I = e,fi, r. 

Neutrinoless double beta decay (Oz//5/5) is known to be a sensitive probe of Majorana 
neutrino masses and mixing. Presently the best the experimental lower bound [15] on the 
0^/?/9-decay half life was obtained for 76 Ge: 

r°; 2 ( 76 Ge) > 1.9 x 10 25 yrs- (39) 

In Ref. (2! this bound was used to constrain the contribution of Majorana neutrinos of 
arbitrary mass. In a good approximation this constrain reads: 

V \ U ^ e%ahm ^ < 5 x 10-8 GeV- 1 . (40) 

with go — 105 MeV. Here denotes the Majorana CP-phase of the Majorana neutrino 
state of mass m u k- For the light-heavy neutrino scenario one has: 

E ^^+g - 2 £ |f/ ei | 2 e^m^<5xlO- 8 GeV- 1 . (41) 
Ar=heavy N «=light 

where 3> go and m ui g . This stringent constraint, applied to each term separately, 
seems to make unrealistic observation of like-sign dielectrons (positrons). However, the 
possibility of CP violation in the neutrino sector should be taken into account. If there are 

10 



m N [GeV] 



FIG. 3: The cross section of p + pip) — > e e jj + X mediated by a heavy Majorana neutrino, at 
the LHC (solid line) and Tevatron (dashed line), as a function of m^. The thick lines correspond 
to the constraint from the Oz^/3/3-decay experimental data and the thin lines are calculated from 
the fixed mixing elements \U t n\ 2 ~ \U^n\ 2 ~ |£4aH 2 = 10~ 3 . 

more than one heavy neutrino state Nk, with different CP phases a at, then different terms in 
the first sum of (jUJ) may compensate each other, reducing the individual Oz//3/3-constrain on 
each of them. Note that even with the presence of only one heavy neutrino there may happen 
the above mentioned compensation between the heavy N and the three light neutrinos 
since they all coherently contribute to (V/3/3-decay. Thus searching for like-sign dielectron 
(positron) events may provide evidence for CP violation in the neutrino sector. In case of 
observation of these events at a rate larger than that stemming from the Oz//3/5-limit applied 
to each term in (I41j) . this would point to CP- violation in the neutrino sector. 

With this in mind we have calculated the maximum of the p + p(p) — > e^e^jjX and 
p + pip) — > e ± e ± W T X cross sections at the Tevatron and LHC, consistent with the Ovftft- 
limit (f4TI) . assuming no CP violation in the neutrino sector. The results are shown by the 
thick lines in Figs [3] and HI In the calculation we used the limits 

< 5 x 10- 8 GeV-\ |LU| 2 <10- 4 , 
\U tN \ 2 < 1(T 4 for 3GeV < m N < 80GeV, (42) 
\U rN \ 2 < 10- 2 for 80GeV > m N . 

For comparison, we also give the cross section calculated with the assumption | fT^-iv 1 2 ~ 
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m N [GeV] 



FIG. 4: Same as in Fig. El but for p + p(p) -► e ± e ± W T X. 



\Ufj,N\ 2 ~ | C/ e AT | 2 = 10" 3 , shown in Figs. [3] and 0] by thin lines. 

With the same constraints (1421) we calculated bounds on p + p(p) — > I I jjX and p + 
Pip) ~ * l^l'^W^X cross sections for other lepton flavors, at the LHC and Tevatron, and the 
results are plotted in Figs. [5] and [61 

In Tables. |T] and [IT] we show the upper limits for p + p(p) — > I I jjX event number 
for all possible lepton flavors, calculated with the constraints in Eq. (|42|) . at the LHC with 
integrated luminosity 10 fb" 1 , and at the Tevatron with integrated luminosity 2fb _1 . 

The numbers given in these tables, except for the like-sign /x/x-production, do not take 
into account event selection cuts necessary to suppress the background and which may also 
dramatically affect theoretical predictions for the signal. In the present paper we do not 
address this issue. In order to evaluate an impact of the experimental cuts on our results 
we use the results of the recent analysis of the like-sign /x/x selection criteria for the 



case of the LHC and TEVATRON experiments. In Ref. [5] there were found sets of the 
cuts which permit efficient suppression of the background and confident selection of /x ± /x ± - 
signal. In Tables I and II we indicate in brackets the number of /x ± /x ± -events after applying 
these cuts for the LHC and Tevatron experiments. As seen from these tables, the regions 
where there is still a chance to observe /x ± /x ± -production with statistical significance are 
20GeV < m N < 70GeV for LHC and 20GeV < m N < 50GeV for Tevatron. As to the other 
lepton flavors there exist in the literature studies (see, e.g. Refs. jl, 8] ) of the corresponding 
selection cuts, although not so detailed as in Ref. Unfortunately the effect of these cuts 
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10 100 10 100 

m N [GeV] m N [GeV] 



FIG. 5: The bounds on the cross section of p + p{p) — ► l^l jjX at the LHC (solid line) and 
Tevatron (dashed line) as a function of rriN for W = tt, [it, e\x and er. 

on our results cannot be directly recognized and requires a special analysis to be done 
elsewhere. Roughly, according to Refs. S], the cuts reduce the like-sign dilepton event 
by an order of magnitude for < lOOGeV, leaving practically unchanged the region of 
higher m N values. 

The following comment is in order here. Notice that because of the involvement of all the 
mixing matrix elements U e N, U^n, U t n in the like-sign dilepton production cross section ( 12(1 - 
( |22|) . one needs to rely on various scenarios relating these matrix elements, in order extract 
limits on them from experimental data. These additional not well motivated assumptions 
decrease reliability of such results. Naturally it would be desirable to get the possibility 
of extracting individual limits on each of the three mixing matrix elements U e j^, U^n, U t n 
without any ad hoc assumptions on their relative sizes. Such an extraction is possible in 
the studied case of production of like-sign dileptons by simultaneous search for events with 
different lepton flavors. As follows from (T20l) - (r22l) 

Thus, searching for one flavor diagonal and two flavor non-diagonal processes and establish- 
ing upper bounds on each cross section within the same kinematical domain, one can set an 
upper limit on the corresponding mixing matrix element \Uw\. Therefore searches for like- 
sign dileptons with all the lepton flavors are required for the model independent extraction of 
each mixing matrix element \Uin\. An experimentally challenging point is searching for the 
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m N [GeV] 



m N [GeV] 



FIG. 6: The bounds on the cross section of p + p(p) -> iH'+W^X at the LHC (left fi ffure) and 
Tevatron (right figure) as a function of tjim for IV = fj,fj,, tt, /xt, e/i and er. 
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TABLE I: Maximum event number of pp — ► I I jjX at the LHC (10 fb _1 data) consistent with 
the limits in Eq. (I42|) . Indicated in brackets are typical number of events after the event selection 
cuts of Refs. [5, 6]. 

events with r's since they involve missing momentum. In this respect it is worth mentioning 
a method for reconstructing r-momenta on the basis of the analysis of the isolated charged 
tracks from r-decays proposed in Ref. 0. The authors have argued their method to be a 
promising basis for future searches for like-sign dileptons with one or two r's. 
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TABLE II: Same as Table. HI but for pp — ► l^l^jjX at the Tevatron (2 fb 1 data). 
V. SUMMARY AND CONCLUSIONS 

We have studied like-sign AL = 2 dilepton inclusive production mediated by heavy Ma- 
jorana neutrinos. We focussed on the dominant mechanism via resonant Majorana neutrino 
production. We applied the existing limits on heavy Majorana neutrino mass and its 
mixing with active flavors u e , z/ M , u T , in order to predict maximal like-sign dilepton produc- 
tion rates ior p + p(p) — > l^l jjX and p + p(p) — > l^l'^W^X , at the Tevatron and the LHC. 
Special emphasis has been made on the stringent limit on the heavy Majorana neutrino 
derived from 0i//3/9-decay experiments. Our results displayed in Figs. [3] -[6] and Tables I, II 
have been obtained assuming no CP-violation in the neutrino sector. 

As seen from Tables I and II, we expect no good prospects for observation of like-sign ee 
events from the studied reactions neither at the Tevatron, with 2 fb -1 integrated luminosity, 
nor at the LHC, with an integrated luminosity of 10 fb -1 . This is a direct consequence of very 
stringent 0^/9/?-decay limits. With an increasing integrated luminosity 100 fb~ x at the LHC, 
the like-sign ee production can reach a few events in the region 20 GeV < < 60 GeV. 
For the p + p(p) — > e ± e ± W T X process, in which a real W boson is detected 100 fb -1 data at 
the LHC are sufficient for discovering Majorana neutrino through like-sign ee production. 
One of the messages of the present paper is that despite of these discouraging predictions 
the like-sign ee events are worth searching for both at Tevatron and LHC. The point is that 
any observation of such events in these experiments would point to the CP-violation in the 
neutrino sector, as explained in sec. [TVl 
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Like-sign dilepton production for other flavors at the LHC and Tevatron have bet- 
ter prospects. In the case of /i ± /i ± production, where we were able to estimate the ef- 
fect of the event selection cuts, the possibility of their observation is open in the regions 
20 GeV < m N < 70 GeV for LHC and 20 GeV < m N < 50 GeV for the Tevatron. 

We proposed a method of model independent extraction of the heavy Majorana neutrino 
mixing Uim with the neutrino active flavors, on the basis of searches of certain sets of like- 
sign dileptons with different lepton flavors. This method will be helpful in reconstructing 
the structure of heavy neutrino sector. 
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Subatomic Physics" and CONICYT No.FB 0821. 



APPENDIX A: INTEGRATION LIMITS FOR s a/b AND t a/b 

The limits of integration over s~, b , s^, b ,t~< b and t£, b , given in Eqs. (I2"0"j) and fT2"2l . can be 
obtained in the following way {{J. The fourth order Gram determinant has the form 



A 4 (-FV, h, k a , k b ) 



Pw 


■ Pw 


Pw 


■h 


Pw 


■ k a 


Pw 


■ h 


k- 


Pw 


k- 


h 


k- 


ka 


k- 


kb 


k a 


Pw 


k a 


h 


k a 


k a 


k a 


h 


h ■ 


Pw 


h ■ 


k 


h ■ 


k a 


h ■ 


h 



(Al) 



From the definitions in Eq. f|T5l) . we see that A 4 is a function of s, s a , s&, t a and t b . Picking 
Sb as the innermost integration variable, explicit evaluation of A4 yields 



I6A4 = asl + bs b + c = a(s b - s b )(s b - s b ). 
Therefore the s b limits are 



(A2) 



s b 



2a 



as h 



(A3) 



The s a -integration limits can be obtained from the requirement A > 0. Solving this equation 
yields 

$t = (h(s + m 2 b - m 2 a ) + s(m 2 b - mfj + m\{m 2 a + m 2 b ) - m 2 b { 



-m b {m b -m a ) 
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s a = t b + m 2 a + ((tl + m 2 - m 2 h ){m 2 l2 -t a - t b ) 

- \ 1 / 2 (t a ,t b ,P^)\ 1 / 2 (t a ,mlml))/(2t a ), (A4) 

where the third order Gram determinant has the form 

X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz. (A5) 

Defining 

Q = s{t b + m 2 b - m 2 h ) - m 2 a t b + m\[m 2 a + m\) - m 2 b {m 2 b + m\ + P^) 

+ ^ 2 (s,m 2 a ,m 2 b )^ /2 (t b ,m 2 b ,m 2 2 ), 
a' = 2m 2 b {Q + m 2 b {t b + m 2 a + m\ + T%)) 

b' = Q 2 - mi(t b - m 2 a -Pl + ml) 2 - 4mt(m 2 a + m 2 h )(t b + P%) 
d = 2m 2 b (Q(m 2 a -m 2 i )(t b -P^) 

+ m 2 b (m 2 a - m 2 h ) 2 {t b + P^) + m\{t h - P^) 2 (m 2 a + mfj), 

we can arrive at the integration limit for t a : 

-b' - \jv 2 - Aa'd 



Finally we give the t^-integration limits as 

4 = (m h + JP 2 ?) 2 + m 2 - ((s + m 2 a - m 2 b )(s + (m h + /p^) 2 - m 2 ) 

± (A 1 / 2 (S,m 2 ,m 6 2 )A 1 / 2 (S,(m il + v ^;) 2 ,m 2 2 )))/(25). (A12) 



(A6) 
(A7) 
(A8) 

(A9) 



(A10) 
(All) 
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